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Scope and Purpose—Planners for spatially distributed service systems must confront the problem 
of locating their service facilities. In this paper we address the situation in which one such facility 
is to be located to service an entire region. The facility is assumed to house a given number of 
mobile service units (e.g. emergency repair vehicles, ambulances), any one of which can be dispatched 
to handle an on-scene service request. Both the arrival times and service requirements (i.e. service 
times) of service requests are probabilistic, representing an advance over earlier models. Queueing 
delays are not allowed, so it is assumed that any service requests that occur while all service units 
are busy are handled (at a cost) by some back-up service system. While our model and analysis 
differ from traditional deterministic location models, we show that the optimum location for the 
facility is the same as that obtained from a well-known deterministic model. 


Abstract—We consider a demand-responsive service system in which n mobile units (servers) are 
garaged at one facility. Service demands arrive in time as a homogenous Poisson process, but are 
located over the service region according to an arbitrary probability law. Given a random service 
demand, either (1) a mobile unit is dispatched to the demand’s location to provide on-scene service 
or (2) the demand is lost (i.e. it is handled by some back-up system). The resultant queueing system 
is an M/G/n loss system operating in steady state. The objective is to locate the garage facility so 
that the average cost of response is minimized, where the cost of response is a weighted sum of 
mean travel time to a random serviced demand and the cost of a lost demand, the weights being 
the respective probabilities of occurrence. We show that the optimum facility location reduces to 
Hakimi’s well-known minisum location. 


1. INTRODUCTION 


The motivation of this paper comes from consideration of facility location problems in 
which spatially dispersed ‘“‘customers” receive service from mobile “servers,” garaged at 
one or more “home locations” while not servicing customers. In the public domain, examples 
include ambulance services and fire departments. In the private domain, important examples 
are emergency repair services and certain home delivery services. 

In many classical locational decisions, minimization of average system-wide travel time 
has been the sole objective. This decision criterion implies that servers are available at all 
times and each call for service by a customer is answered immediately by the nearest (always 
available) server. Under such a congestion free environment, average travel time is a good 
measure of system performance. In an urban environment, service-oriented systems are 
plagued with uncertainties: temporal and spatial uncertainties as to the arrivel time and 
location, respectively, of the next customer; service time uncertainties as to the amount of 
time required to complete a service; travel time uncertainties as to the fluctuations of transit 
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time. As a consequence of such uncertainties, in urban emergency service settings where 
system utilization (i.e. the fraction of time that a server is busy) is high, customers often 
find all servers unavailable at the time a service request is initiated. Such a customer could 
wait in a queue until an appropriate server is free, or the request for service could be lost 
at a cost (e.g. being served by a distant or more expensive back-up unit). It is the latter 
case, the case of “‘loss,”’ that we will analyze in this paper. 

Most previous queueing considerations in this context relied on the assumption of 
exponentially distributed service times. Such an assumption allows one to model the system 
as a continuous time Markov process. Mirchandani, Silva, and Visocki [10] study the 
resultant Markov model for a network with two nodes. Berman and Larson [1] extend the 
Markov analysis to an arbitrary network having servers, any combination of which may 
be unavailable (i.e. busy) at a given moment. They obtain Hakimi-type [5] nodal optimality 
results; i.e. there always exists a set of optimal locations on the nodes of a network. They 
assume that the state probabilities of server status (free or busy) remain fixed for local 
changes of service facility location (except when the change of facility locations results in 
changes of server assignment policies). This approximation is good only when on-scene 
service time is much larger than travel time. 

In many applications of interest, the total service time (defined to be the sum of travel 
time and on-scene [and perhaps related off-scene] time), is dominated by travel time. The 
resultant distribution for total service time is (i) non-Markovian [i.e. not exponential] and 
(ii) dependent on the location(s) of the server(s). Allowing for such complications, Berman, 
Larson, and Chiu [2] prove the nodal optimality result for the one server loss case (demand 
is lost at a cost when the lone server is busy). They do not make any restrictive assumption 
about service time distribution nor do they assume constancy of the state probabilities upon 
small changes of facility locations. 

In this paper, we extend the model of Berman, Larson, and Chiu [2] in two ways: (i) 
we wish to locate a single facility for » servers; and (ii) we allow any topological structure, 
distance metric, and demand distribution in the loss model. When our result is applied to 
a network structure with only nodal demands, it is both more general and more restrictive 
than the problem studied in Berman and Larson [1]. It is more general because we make 
no assumptions on service time distribution and we allow the steady state probabilities to 
vary continuously as we change facility location. It is more restrictive because we are 
locating one single home for all servers. 

Our objective is to locate a single j..age facility for n servers, which will minimize 
the expected weighted cost of travel time and the cost of lost customers. We will call such 
a location an n-server-single-facility-loss median (n-SFLM). We will show that the standard 
average travel time minimizing location coincides with our n-SFLM under any topological 
setting with any demand distribution over the region of interest. The situations include 
continuous link demands on a general network (see Chiu [4]), location problems on an 
Euclidean plane, on a plane with L, (or L,) metric (see Larson and Sadiq [8], and Odoni 
and Sadiq [11]). The cost of “loss” is, curiously enough, only constrained to be nonnegative. 
Section 3 contains a discussion of the behavior of our objective function when the underlying 
topology is a network with discrete nodal demands. We end this paper with a discussion 
of the difficulties one encounters in generalizing the loss model we study. 

For general background on facility location problems, we refer the reader to Francis 
and White [6], Thisse and Zoller [13], and Francis and Mirchandani [9]. 


2. MODEL FORMULATION AND RESULTS 


Consider a region R in which customers’ requests for service are generated as a time 
homogeneous Poisson process with rate A. The location of such a service request is governed 
by some probability distribution function over the region of interest. A single server, if at 
least one of the total pool of nv is available, will travel from the home location of the n 
servers to the point of service request, provide the required service, return to the home 
location, and prepare itself for another request. Figure 1 shows the temporal sequence 
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d(x,y) = Travel distance between x and y 
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Fig. 1. Temporal sequence of a service duration for a facility at x and incident at y. 


comprising a single total service time, during which the server is considered busy. A request 
for service is lost at a nonnegative cost Q if all ” servers are busy upon its arrival. Our 
objective is to locate the garage facility for the n servers so that the weighted sum of mean 
travel time (when at least one server is idle) and the cost of loss Q (when all servers are 
busy) is minimized. 

Before we analyze and search for the n-SFLM, we make a few remarks: (1) For a 
garage facility located at x, the service time distribution depends only on call location and 
not on system status, nor does it depend on server identity. (2) We are locating one facility 
for n servers. (3) The location of the facility determines not only the service time distribution 
for each customer location, it also determines the system’s steady state probabilities. (4) 
For a facility located at x, the objective function takes on the following form: 


Z(x) = [1 — PCOod]tx) + PDQ, (1) 
where 


P(x) = probability of saturation (i.e. all servers busy), 
1(x) = expected one-way travel time, 


Q = cost per lost customer, Q 2 0. 


The argument x denotes location of the facility in the region of interest. The expected travel 
time can be viewed as expected response time, since a customer is either served immediately 
or lost. (5) All servers are indistinguishable in terms of their service time distribution and 
hence the expected service time. To evaluate Z(x), we need to know the expression for 
P(x). This is readily available as the Erlang Loss Formula (Takacs [12]): 


P{x) = a eA (2) 
Dolev 
p(x) = Aw + ABt(x) = ASs(x), (3) 


s(x) = w+ Br(x) 
= expected service time for one service request, 


expected nontravel related component of total service time, 
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B => 1 accounts for round trip travel time; e.g. when B = 
it takes as much time traveling from the facility at x to the 
incident location as it does on the return trip. 


Thus $(x) consists of round trip travel time and other nontravel related service time. 
Our optimization problem becomes 


Min Z(x), 
xeR 


where R is the region of interest. 
We define a point x in R that minimizes ¢(x) as a minisum location. We will first 
prove a weaker version of our main result equating the minisum location with the n-SFLM. 


THEOREM I: 
The n-SFLM coincides with the minisum location if Q@ > t(x) for all x in R. 
Proof. Since the minisum location minimizes ¢(x), we need only show that 


ae >= 0 for all xeR, 
but 
G20) 2 PGs aRO = )] oF oe > 0 
d(x) 
if 
dP.) > 0, since [1 — P,(x)], A, B, and [Q — t(x)] are all nonnegative. 
dp(x) 
since 


tee aby (x)ili iF: p(x)" apy (x (fr ‘) > 0, our result follows. 


dp(x — 1)! 


We offer the following intuitive interpretation for Theorem 1: to minimize Z(x), we 
have to minimize P,(x), the probability of saturation. This is because the cost of loss Q is 
larger than the travel time #(x). Thus we would like to minimize P,(x), the chance of system 
saturation. P,(x) is the probability that all servers are busy and p(x) is the utilization factor 
which is linear in #(x) (since W is independent of server location). Increasing ¢(x) increases 
the expected service time for each call. This, in turn, should produce more work for the 
servers and thus increase the system saturation probability. Therefore, a location that 
minimizes ¢(x) also minimizes Z(x). 

It turns out that we can relax the bound on Q and obtain the same result. The proof 
is direct but algebraically involved. 


THEOREM 2: 
The n-SFLM coincides with the minisum location for Q > 0. Proof: For simplicity, 
we will suppress the argument x in P,(x), r(x), etc. 


eG oF) tee > (1 — P) 


dP, 
dt a : dp p dp 


> 


since ABQ(dP,/dp) > O from Theorem 1, and p = AW + Afi = At. 
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Therefore, we need only to show that (1 — P.) — p(dP,/dp) = 0 


dP, 
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After some algebraic manipulation, B is seen to be 


Since each term in A is positive, we need only compare terms in A that involve p"+* with 
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We wish to bound the denominator of each term of the above expression by n!k!. That 
is, we wish to show that n!k! > (mn — )!(k + J! fori < n — k. This is true if and 
only if 


nin —1)...2—74+ Dn — DIKE > mM —DUk+ Dk+i-1...(k+ Dk! 


or 


n(n — I(n — 2)...@ —14+ 1 > Kk + DA+i-1)...(k + 1) after cancelling 
similar terms on both sides of the inequality. 


Comparing the above inequality term by term, ie. (n — J) vs (kK + i- JD, / = 0,1, 
2,...,— 1,wehaven—/>k+i—I/sinceen—k >i. 
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Thus, 


dP as 
dp 


(1 — P) — Apr 
and 
dZ 
— 2 0 follows. 
dt 
Therefore a location minimizing t(x) also minimizes Z(x). 


3. SPECIALIZING TO A NETWORK TOPOLOGY 


In this section, we will discuss the behavior of Z(x) on a link of a general network 
with discrete nodal demands. Traveling between any two points on the network takes on 
the shortest path. In order to minimize the amount of new notation, we will state our 
results without proof. We refer the interested readers to Chiu [3]. 

It is well known that f(x) is piecewise linear and concave on a link of a general network 
(e.g. see Hakimi [5]), where x denotes location of the facility on the link. This results in 
the nodal optimality condition for a point minimizing t(x). We will define the region over 
which ¢(x) is continuously differentiable on a link as a primary region, and the boundary 
between two adjacent primary regions as a breakpoint. Figure 2 shows the behavior of t(x) 
on link (ij) of a general network. 

After explicitly expressing f(x) in terms of shortest path travel distances between 
demand nodes and facility at x, one can easily prove the following results (see Chiu [3] for 
details): 


LEMMA | 

Z(x) is monotone over any primary region on a link of a general network. 

Since Z(x) is nondifferentiable across breakpoints [as in the case of f(x)], one has to 
evaluate the directional derivatives of Z(x) at a breakpoint from both directions in order 
to examine the behavior of Z(x) on a link. To characterize the behavior of Z(x) across a 
breakpoint one has the following result: 


LEMMA 2 


The directional derivative of Z(x) at a breakpoint is nonincreasing. Lemmas | and 2 
are sufficient to guarantee that at least one minimizer of Z(x) is at a nodal location. 


t(x) 


Primary << t x —> 
region breakpoint 


@) 


Fig. 2. ¢(x), primary region and breakpoint on link (i, /). 
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Furthermore, if we assume that Q > ¢(x) for all x on the network, one can prove the 
following: 


LEMMA 3 

If Q > t(x) for all x on the network, then, (i) Z(x) is monotone and concave on any 
primary region and (ii) Z(x) is concave on any link of the network. 

However, Lemma 3 is not required to prove a Hakimi-type nodal optimality result. 
These three lemmas are quite immediate given Theorems 1 and 2, and the behavior of t(x) 
on a link. We will close this paper with a general discussion of the difficulties one encounters 
in generalizing the loss model we just analyzed. 


4. CONCLUSION AND DISCUSSION 


We have formulated the n-server-single-facility-loss-median problem in which we want 
to locate a single facility to house n servers in order to minimize the weighted average of 
mean travel time and the cost of lost customers. Servers are indistinguishable as far as their 
service time distributions are concerned. This allows us to use the Erlang Loss Formula. 
We do not have to evaluate the steady state probability of each combination of server status 
(busy or idle). Only the knowledge of the saturation probability is required. Otherwise, we 
have to resort to numerical methods to compute each state probability as discussed in Wolff 
and Wrighton [14]. 

We are able to prove the equivalence of the minisum location and the n-SFLM without 
imposing any lower bound on the value of Q (except for nonnegativity). This is very peculiar 
since the cost for a lost customer can take on the value zero. Intuitively, to minimize Z(x), 
one has to strike a balance between 7¢(x) and P,(x). When we place the facility at a very 
“bad” location, the expected service time increases due to the increase of travel time 7(x). 
This will increase the saturation probability P.Qc). Even if Q is zero, in which case there is 
no penalty for losing a customer, the increase in P(x) [decreasing the chance that some 
server(s) is available] is countered by the deterioration of the travel time [and hence service 
time s(x)]. The final mathematical analysis shows that the best location is still at a place 
where ¢(x) is minimized. 

Within the framework of our model and the subsequent analysis, one can assign a 
server preference list (i.e. dispatch policy) at each demand point without violating any of 
the conclusions. Such a list allows preferences for bilingual personnel or servers familiar 
with local neighborhoods. Such generality also ties in with previous work of Larson [7]. 

We have contemplated several extensions to this loss model. The first obvious extension 
is that servers are allowed to reside at different locations. This modification gives rise to 
different service time distribution for different servers. Hence, the Erlang Loss Formula 
will no longer be valid. A second extension is to allow different types of calls to be handled 
differently by different servers and thus model w, the nontravel component of total service 
time, to be server specific. This, of course, leads to the same difficulty. Another modification 
is to allow a server to respond to another service request as soon as it completes on-scene 
service. This modification leads to statistical dependence between successive service times, 
and the difficulty of deciding a dispatch policy. Available queueing theory results cannot 
handle the above difficulties. 
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